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We model the charge, spin, and heat currents in ferromagnetic metal|normal metal|normal metal 
trilayer structures in the two current model, taking into account bulk and interface thermoelectric 
properties as well as Joule heating. Results include the temperature distribution as well as resistance- 
current curves that reproduce the observed shifted parabolic characteristics. Thin tunneling barriers 
can enhance the apparent Peltier cooling. The model agrees with experimental results for wide 
multilayer pillars, but the giant effects observed for diameters < 100 nm are still under discussion. 


I. INTRODUCTION 

Thermoelectric effects result from the coupling be¬ 
tween energy and particle transport in conductors. An 
example is the heat current that is associated with a 
charge current and proportional to a material parameters 
called Peltier constant. At a thermocouple, i.e. a junc¬ 
tion between two conductors with different Peltier coef¬ 
ficients, the heat current is not conserved, which implies 
heating or cooling depending on the current direction.b^ 
The thermopower, on the other hand, is the thermoelec¬ 
tric voltage that is generated by a temperature difference 
over a metal wire that is proportional to the Seebeck co¬ 
efficient. A thermocouple generates an isothermal ther¬ 
moelectric voltage proportional to the difference between 
the Seebeck coefficients when the end of the wires are 
maintained at a temperature different from the junction. 
Nanostructured materials can enhance the efficiency of 
thermoelectric devices;^ Thermoelectric effects in metal¬ 
lic heterostructures including ferromagnets depend on 
the spin degree of freedom^ The spin dependence of 
thermoelectric cooling^^— is part of the field that studies 
the coupling between spin, heat and electric transport 
in small structures and devices, or spin caloritronicsi^ 
Heating is an important issue for spin torque magnetic 
random access memories (STT-MRAM) device [MRAM], 
and spin caloritronic effects can improve their perfor¬ 
mance. 

An enhanced Peltier effect has been reported by 
Fukushima et alr^^— in metallic multilayers when struc¬ 
tured into nanopillars. The effect was detected by change 
AR of the resistance Rq as a function of current bias 
that acted as a thermometer. The Peltier effect cools or 
heats the systems by a term linear to the applied charge 
current Ic and Peltier coefficient H, while the Joule 
heating induces a temperature and resistance change 
that scales like the square of the current bias such that 
AR ^ RqIc ~ nJc- At small currents the linear term 
dominates and causes a reduction of the resistance, i.e. 


an effective cooling, that in some structures was found to 
be very large. The Peltier coefficient was found by mea¬ 
suring the current where heating and cooling compensate 
each other and AR = 0 and therefore H = RqIc^K 

The observed H = 480 mV in pillars containing Constan- 
tan is attractive for cooling nanoelectronic devices.— The 
cooling power enhancement was tentatively explained 
by Yoshida et a/— by adiabatic spin-entropy expan¬ 
sion. However, such an equilibrium cooling mechanism 
could not explain that H is material dependent and even 
changes sign. The diffusion equation approach by Hatami 
et ah^ did take not into account either the precise sample 
configuration nor Joule heating and could not reproduce 
the large observed effects. The physical mechanism of the 
giant Peltier effect therefore remains unexplained. On 
the other hand, the recent experiments by Bosu et un¬ 
confirmed large Peltier coefficients for pillars including 
Heusler alloys when becoming very narrow. The present 
research has been motivated by the wish to model the 
heat and charge currents realistically in the hope to shed 
light onto this quandary. We report detailed calculations 
for the structure and model parameters matching Bosu et 
experiments and compare results of semi-analytic 
calculations with experiments. This study is limited to 
thermoelectric effects as described by the two-current 
model of thermoelectric transport in which spin current 
is carried by particle currents. We do not include explic¬ 
itly phonon contributions to the heat current as well as 
phonon/magnon drag effects on the thermoelectric coef¬ 
ficients, which may lead to a temperature dependence of 
the model parameters. Furthermore, we completely dis¬ 
regarding collective effects that give rise to e.g. the spin 
Seebeck and spin Peltier effectsi^ There are no indications 
that these approximations will do more than leading to 
some renormalization of the model parameters. While we 
are still far off a complete understanding of the experi¬ 
ments, we find evidence that very thin (Ohmic) tunnel 
junctions can enhance the Peltier effect. 

This paper is organized as follows. In Section H, we re- 
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view the standard Valet-Fert model for spin transport^ 
in our nanopillars, with explicit to inclusion of interfaces. 
In Section III, we extend the model to include heat cur¬ 
rents, charge and spin Joule heating, and explain our 
method to compute temperature prohles. In Section IV, 
we present results for the Peltier effect due different inter¬ 
facial thermoelectric parameters and simulations of the 
Peltier effect are also performed, illustrating the impor¬ 
tance of interface resistances, to finish in section V with 
a summary and conclusions. 


II. SPIN-DEPENDENT DIFFUSION IN F|N|Nb 
MODEL 


strong, fi approaches the Fermi-Dirac distribution which 
depends on temperatures Ti and chemical potentials jii. 
We disregard spin-dependent temperatures here^ but al¬ 
low for spin accumulations, i.e. local differences between 
chemical potentials for both spins. 

The spin particle and heat currents at a po¬ 
sition cc in a resistive element are^S 

=A,|dej(“)(e,x) (1) 

and 

dee/“)(e,x)-/xo J (2) 


Our model can be applied quite generally to arbitrary 
multilayered structures, but we focus here on the charge- 
current biased trilayer nanostructures measured by Bosu 
et alr^ that are composed of a ferromagnetic metal F and 
two normal metals N and Nb, respectively, as sketched 
in Fig. [I] The thicknesses of F, N and Nb are Lp, L 
and Lb, consecutively, and the device is sandwiched be¬ 
tween two thermal reservoirs at same temperature Tq. 
The electric, spin and heat transport is described by 
an extended Valet-Fert model^ including interfaces'^ 
and spin-dependent thermoelectric effectsi^ The param¬ 
eters are interfaces resistances i?i and i ?2 for interfaces 
the F|N and N|Ng respectively ; bulk resistance Ri 
[i = F,N, Nb) for each metal, as well as the spin polar¬ 
ization Pp of the ferromagnetic metal. 

We adopt a one-dimensional diffusion model in which 
the currents flow along the x-direction and the origin is 
at the F|N interface. In the collinear two-channel resistor 
model, the electrons are in either spin-up or spin-down 
states. We divide the structure into various elements such 
as resistors, nodes and reservoirs. Discrete resistive ele¬ 
ments are interfaces, tunnel barriers or constrictions that 
limit the transport. For our purpose, resistors are sepa¬ 
rated by nodes in which electrons can be described semi- 
classically by distribution functions fi. If the interac¬ 
tions electron-electron or electron-phonon are sufficiently 
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FIG. 1. (Color online) F|N|Ng nanopillar structure biased 
with a charge current R. We consider a nanopillar with cross 
section Ac- F, N and Nb represent a ferromagnetic metal, 
normal metal spacer and normal metal B respectively. Lp, L 
and Lb are the thicknesses of each section. The ends of the 
nanostructure are connected to thermal reservoirs kept at a 
constant temperature Tq. 


respectively, where is the cross sectional area of the 
nanopillar, a =t (4.) is the electron spin degree of free¬ 
dom, is the spin, energy (e), and position (x) depen¬ 
dent spectral current density, and /xq the ground-state 
chemical potential. (e, x) is described 

by local Fermi-Dirac distributions f^°‘' at temperature 
T and spin-dependent chemical potentials times the 
energy-dependent conductivity a°‘{e). 

The spin accumulation is defined as fig = Ft ~ 
where are the spin-up and spin-down chemical po¬ 
tential of the material, while the charge chemical poten¬ 
tial is the average of the sum of spin-up and spin-down 
chemical potentials /Xc = {^t + Mi)/2- 

The transport in each layer is governed by spin and 
charge diffusion equations^ given by 


9^ _ Fs 

dx^^" ^^2A2 


( 4 ) 


where A is the spin-flip diffusion length, usually much 
larger in normal metals than in ferromagnetic metals 
^n,Nb Af, and 


Pf 


t I 
o-p - o-p 

t , i 

(7p J- (7p 


( 5 ) 


which is the spin polarization of conductivity in terms of 
the spin-dependent conductivity for each channel in the 
ferromagnet. In normal metals N and Nb these polariza¬ 
tions vanish {Pn,Nb = 0). The solution of Eq. ([3]) 


fj-s = Ye^ + Ze"^ 


( 6 ) 


depends on the integration constants Y and Z. For a 
ferromagnetic metal F,— we find (see Fig. [IJ 


F^f\x) 


Ic 

apAc 


O’ fP X 

x + ai± —Ae^F 

CTp 


( 7 ) 


where ai is the voltage drop at the interface F|N, up 
is the electrical conductivity of the ferromagnetic metal. 
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cr^ the spin-dependent conductivity, Xp the spin-flip dif¬ 
fusion length and /3i a coefficient to be determined by 
boundary conditions at the interface.— 

For the normal metals, the spin-up and spin-down 
chemical potentials read 

^ =-^ 2 ; ± bie ± b2e^N ( 8 ) 

e ctnAc 


Material 

A (nm) 

p (/iflcm) 

C 02 MnSi 

2.1 

70.0 

Au 

60 

2.27 

Cu 

350 

1.73 


TABLE I. Spin-flip diffusion length and electrical resistivity 
at 300K used for the F|N|Ng nanopillar structure . 


Ic f 

-=- —[x — L)+a2±cie 

e aNsAc 

X — L 

±C2e^ ( 9 ) 

for N and Nb respectively, where 61 , 62 , ci and C 2 com¬ 
plete the number of coefficients that describe the spin- 
dependent transport in the present trilayer system. The 
spin accumulation in each layer of the F|N|Ng nanowire 
are while the charge chemical potentials read 

fj,^{x), and the spin-dependent current^^ in a bulk ferro¬ 
magnetic metal is (Ohm’s Law): 


B. Boundary conditions 

The boundary conditions are spin and charge current 
conservation at the interfaces. 

rAo) = r^{o) = ( 13 ) 

for the F|N interface and 

lAL) = IA{L) = I 2 (14) 

for the NjNg interface. We assume that the spin accu¬ 
mulation vanishes at the end of Nb 


Ix\x) = -Aca 






( 10 ) 


where X = F,N, Nb and = (TAr/2. The spin current 
(t) (-i) 

= /jj. — is the difference between spin-up and 
spin-down currents 

where parameter such as = pxXxjAf.^ which is 
the resistance over the spin-flip diffusion length Xx in 
X and px is the corresponding electrical resistivity, are 
implicit in the calculations. 


— 0 ( 15 ) 

which is valid for Xx or Lb ^ Xnb and/or when 

nanopillar diameter widens at Lb- We can now deter¬ 
mine /3i, 6 i, 62 , Cl and C 2 in terms of the coefficients. 

Then, it can be now computed the spin accumulation, 
spin current and charge chemical potential. 

The total electrical resistance, R = pc/ {sic), of the 
device can now be written as 

R = Rf\n\Nb = Rpix = —Lf)—Rxb{x = L+Lb) (16) 


A. Interface resistances 

Next we consider spin-dependent transport through 
the interfaces. We disregard interface-induced spin- 
flips so at the F|N interface— 


RiTo) = 


2Pf 

Ic{l - Pi) 





PfLf 


02 PNbRb 

Ic A, 


(17) 


where 


r(a) 

-'1 



[Li“^( 0 )-M!v“^( 0 )] 


( 11 ) 


oi — IcRi 


2A(Pi-Pf) 

(1 - R^) 


-h Pi{bi -f 62) 


(18) 


where is the interface conductance with polarization 
Pi = /Gi and Gi = G^l^ -b G^. At the 

interface between the two normal metals N|Ng 

/(“) = ^[^W(p)-^M(p)] ( 12 ) 

Charge = /i ,2 = 1^2 + 41 ) 2 ) ^'^‘1 

2 = Ii ^2 ~ ^ 1 ^ 2 ) currents are conserved in the 
interfaces 1 and 2, and assuming that Pi = 1/Gi and 
P 2 = 1 / G2- 


and 

02 = -IcR2 - (19) 

■^C 

are the voltage drop at the two interfaces. 

Numerical results for the transport properties require 
the parameters of the samples considered by Bosu et aZ— 
at room temperature Tq. The ferromagnetic metal is typ¬ 
ically a Heusler alloy Co 2 MnSi (CMS))^ while the nor¬ 
mal metal N is gold and Nb is Cu. The resistivities and 
spin-flip diffusion lengths are given in Table HI 

Fig. [5] illustrates that a charge current R leads to a 
spin accumulation over the spin-flip diffusion length Xf 
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FIG. 2. (Color online) Spin accumulation in a 
Co2MnSi(CMS)[40nm]|Au[10nm]|Cu[100nm] nanopillar with 
diameter D = 180 nm for a current 7c = 5 mA and tak¬ 
ing interface resistances^*^ into account AcRi = 0.915 x 
10“^®flm^,— and AcRi = 3.40 x 10“^®flm^ at 300K, the po¬ 
larization of ferromagnetic metal is Pf = 0.71 and polariza¬ 
tions of the interfaces^ are Pi = 0.77 and P 2 = 0. The dashed 
line shows the spin accumulation when interface resistances 
AcRi and AcR 2 are set to zero (metallic contact). 


/,[mA] 



FIG. 3. (Color online) Spin current in the 

CMS[40nm]|Au[10nm]|Cu[100nm] nanopillar structure 
for a charge current bias 7c = 5mA and taking interface 
resistances into account. The dashed line is the spin current 
when interface resistances A^Ri and AcR 2 are set to zero 
(metallic contact). 


in F, reaching its maximum value at the F|N interface, 
where the spin is injected^ '^^i^^'^° and decays exponen¬ 
tially along the spin-flip diffusion length of the normal 
metals Xn.Nb- The spin current is plotted in Fig. [3] In a 
normal metal the spin current is proportional to the gra¬ 
dient of the spin accumulation, Eq. (fTUl) . It is observed 
in the model (Fig. [3]) that the spin current decays rapidly 
in the central island N. Its behavior depends strongly on 
the spin-flip diffusion length of the metal, for our model 
we have Xau < Xcu- Additionally, it has an influence 
from the design length L,Lb and the boundary condi¬ 
tion established in Eq. m- 


III. SPIN-DEPENDENT 
THERMOELECTRICITY OF F|N|Nj3 PILLARS 


In the experiments the electrical resistance change is 
measured as a function of applied current, reflecting the 
balance between the Joule heating and Peltier cooling. In 
order to model this effect we need to compute the tem¬ 
perature profile distribution T{x) over FjNlNg pillars. 
Temperature distributions have been previously calcu¬ 
lated, but without taking Joule heating into account in 
spin-dependent systemsj^ Assuming that we know the 
temperature dependence of the electrical resistivity p(T) 
and interface resistances i?i _2 (T), the total temperature 
dependent resistance reads 

AR=jJ R [T{x)] dx - R{To) (20) 

where 77(To) is given in Eq. (1171) . For simplicity, we disre¬ 
gard the heat leaked through the cladding of the nanopil¬ 
lar, which is valid when the thermal contact is weak or 
the cladding material has a much smaller heat conduc¬ 
tivity. Significant heat leakage would reduces the tem¬ 
perature gradients calculated here, leading to an overes¬ 
timate of the thermoelectric cooling power. In the follow¬ 
ing we determine the heat current and its divergence in 
the nanopillar taking into account the Kapitza thermal 
resistances at interfaces'^ The temperature profile dis¬ 
tribution along the nanopillar structure is calculated us¬ 
ing heat conservation at interfaces, to finally describe the 
performance of the nanodevice in the resistance-current 
{R-I) characteristics. Except for the temperature depen¬ 
dence of the resistance that serves as a thermometer, we 
disregard the (for elemental metals) weak temperature 
and voltage dependences of the thermoelectric parame¬ 
ters. 

In the Sommerfeld approximation the linear response 
relations between currents and forces in bulk materials 
readi^ 


f Jc\ (I Pf ST \ ( \ 

J« =(T Pf 1 P'pST -d^vLsjile) 

\JgJ \ST P'pST KT/a ) \ -9,lnT ) 

( 21 ) 

where S is the (charge) Seebeck coefficient, a the elec¬ 
trical conductivity, k the thermal conductivity, all at the 
Fermi energy and T is the temperature (disregarding spin 
temperatures^). Here, Jc = /cMc etc., are current den¬ 
sities. 



/f (4 4) 


( 22 ) 


is the spin polarization of the energy derivative of the 
conductivity at the Fermi energy, which is related to the 
spin polarization of the thermopower as 


p _ S\ — SI _ Pp — Pf 
S’-!- -I- 5^ 1 + PpPp 


( 23 ) 
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Joule heating is a source term that causes a divergence 
in the heat current)^ 


_^T - - T 


(24) 


Including the dissipation due to spin relaxation^i^ we 
obtain the matrix expression for the divergence of the 
current densities 



(25) 


T rp rpf rp rp f rp 

L ‘N1 IN1 ‘N2 ‘N2 ‘R 
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FIG. 4. (Color online) Definition of temperatures and 
heat currents in F’IA^IA^b nanopillar structures along the x- 
direction, including Kapitza thermal conductances Gthi, 2 - 


becomes 


A. Heat currents and temperature profiles in the 
bulk of the layers 


The divergence of the heat current in the ferromagnet 
F reads (Eq. (I^ i 


d .2 (I--Pf)m? 

— J, -J,PF + (26) 


which equals the derivative of the heat current in Eq. 

m 


d JF d {Pp- Pf)SfTps d ^ 


{Pp - Pf)SfT Ps 

PF 2 A 2 '"^9x2 ’ 


(27) 

(28) 




kfPf(1 - Pp) 

/ 3 ? 23: 

Pi [e^ - 1 ] + Tni 


2kfA'^ 


kfPf (1 - Pi) 

^ frr. rr. , 2 {P'p-Pf)SfTP^,^^ 


Pi 


2kfAI kfPf (1 - Pi) 


[e^ - 1 ] 


(30) 


and 


d 


^ 2{Pp - Pf)SfTP, , [e^ - 1] 

= " rfPf{1-PI) I" 


Pi 


2e^F [e - 1] 


ffpf (1 — Pi) \ Af 


Lf 


PpF f Lf\ ^ Tl-Tmi 

7F^X-—)*\—Fr- 


( 31 ) 


leading to the heat diffusion equation 


Einally, the heat current distribution reads: 


J. {p'f - Pf)SfT Ps JIpI (i - PDpI 
9a;2 pn 2A2 pn ApnX^ ' 

(29) 

Heat transport is carried in parallel by phonons and 
electrons)^ We assume here efficient thermalization in 
and between both subsystems, meaning that the elec¬ 
tron and phonon temperatures are taken to be identical. 
The total thermal conductivity then reads k = /Cg -f Kp. 

For the ferromagnetic metal F we set T{x = 0) = Tjvi 
where Tyi is governed by the boundary condition at 
the F|N interface discussed in the next subsection, while 
T{x = —Lp) = Tp is fixed by the reservoir (see Fig. |4]). 
The solntion of the heat diffusion Eq. (0^1) disregard¬ 
ing the generilized Thomson effect term then 


• Ag = SfTI, 


2iPp - Pf)SfTPi M - e^ 


(1 - Pi) 


Rf 


Pi 


2e^F 1 _ e 


(1 — Pi) \ R\p Rf 


X Rf\ kfAc 


{Tl-Tni), 


(32) 


where Rp = ppLp/Ac is the electrical and R\p = 
PfXf/Ac the spin resistance. 

Repeating this analysis for normal metals, we obtain a 
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heat current in N 


Substituting Eqs. (1321) and (l33t for a; = 0 leads to 



and Nb 


(x) ■ A, = SnsTIc + (^x - 


‘IL -\- Lb 


- Cl e 


2 i?A. 


- 1 ) 


^Rnb 


2 ^ 2 ' ^ 
e e ^’^B — 1 


2 -Ranb '^Rnb 

C 1 C 2 /2a: 2L + 

R^Nn \^Nb ^Nb 


^Nb 

Lb 


(Tr — Tpf2^ 

(34) 


+ 


KnAc 

L 
2 

i? 2 (l-P|) 1 RTb 


= { - [SnIc 

PI 


KfAcTb IcRf 
' H 2 LF 


2 H 2 


^L, p 

e — 1 
Rf 


fn D ^ F,fAcTf Kjs[AcTfi2 

-^{Rf- Rn) - —^ ^- 




2L fp 

e — 1 


(1 - P 2 ) 


bl 


1 


Rf 

- 1 ]^ 


-h 


j^ 2 i?Ajv 

2R\^ 


+ 


4i?Ar j 

2L \ 

4Rn J i?Ajv \^N 

r, r , '^nAc\ f SfR , KfA. 

S«‘-: + —)[^ + JFL, 

2{Pp - Pf)SfPi 
^ 2 ( 1 -Pl) 

kfAc 



— SfIc — J 
Lf 

2{P'p - Pf)SfPi 

(1 - Pi) 


[e^ - 1 ] 

Rf 


(36) 


B. Interfaces 


Finally, we knit the solutions for the bulk layers to¬ 
gether at the interfaces by boundary conditions. The con¬ 
tacts to an abruptly widening nanopillar may be treated 
as ideal reservoirs (heat and spin sinks) at constant tem¬ 
peratures Tl = Tr = To (see Fig. m. By disregard¬ 
ing interface-induced spin-flipsi^ and, for the moment, 
the Joule heating by the interface resistance, we may 
impose charge, spin and energy conservation at each 
interface /^/°/^ such as Jg {x = 0) = Jg = Jg {x = 0) 
for F|N, where analogous to Fq. (|2T]) 4^ 


and 


SfR kfAc 2{Pp — Pf)Sf/3i 
H 2 LF 



H 2 LF 


2H, 


H2(1 - Pi) 


R\f 


[e - 1] 
Rf 


(37) 


, aJP 

JPA, = GthiA^AT -GiSiTiAfj,^l'> - PpGiSiTi-^ 

(35) 

is the interface heat current, C?thi the Kapitza thermal 
conductance (including the phonon contribution), Ac the 
cross sectional area of the nanopillar, AT = T^i — T'ni 
the temperature drop over the interface, Ti = (T/vi -I- 
Tj(,^)/2 the interface temperature, Gi the electrical in¬ 
terface conductance, the interface thermopower, and 
AMc(s) charge (spin) accumulation differences over 
the interface. 


where Hi( 2 ) = —GiSiA^^P/2 — P'pGiSiA^^P/A ± 
Gthi^c- We may determine the temperatures Tfi 2 and 
T ^2 interface NjNg analogously. 

Fqs. (l36l) and (l37l) include bulk and interfacial Peltier 
effects as well as Joule heating in the bulk materials (see 
Fig.H but not yet the interfacial Joule heating. Here we 
focus on Joule heating by the NjNg interface, which is the 
dirty one in existing experiments. We can treat interface 
heating easily in two limiting cases. In the dirty limit 
the interface is a resistor with small but finite thickness 
Li around the position x = di in which the electrons 
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dissipate their energy directly to the lattice: 

A 7/ = / fr for <^-di< Li/2 

dx \ 0 otherwise. 


(38) 


Clean interfaces, point contacts or coherent tunnel junc¬ 
tions, on the other hand, inject hot electrons (and holes) 
into the neighboring layers where they loose their excess 
energy on the scale of the electron-phonon thermalization 
length In normal metals like Cu it is surprisingly 
large even at room temperature, i.e. A^^ = 60 nmi2^ In 
the clean limit (assuming that A(if -|- A^ is smaller than 
the pillar length) 


dx 


t 2 RrA 
0 


for 


—Ajif < X — di < Ag 

otherwise 


(39) 


The two limits therefore differ only by the volume in 
which the heat is produced. In the extreme case of S> 
Lx all interface Joule heating occurs in the reservoirs, 
where its effect can be disregarded. In the following we 
consider both extremes, i.e. the dissipation occurs either 
in the interfacial thickness L/ or in the reservoirs + 
Af = 00. 

We can implement these models into Eqs. (1551) and 
dSZl) as follows. In Eqs. (1551) and (1551) . Joule heating 
is represented by the power density J^RjAjLi in the 
volume V — ALj. The total power dissipated at the in¬ 
terface is therefore I^Ri- This term can be added to Eq. 
dMD; the hrst term of the third line expresses the balance 
between the Joule heating of the bulk metals to which 
the interface contribution may be added. The interfacial 
Joule heating thereby reduces the cooling power of the 
nanopillar. By contrast, in the ballistic limit and long 
relaxation lengths Joule heating is deferred to the heat 
sinks, and does not contribute at all. In Eq. (|37p the in¬ 
terfacial Joule heating is indirectly related by the already 
determined term Txi of Eq. (1551) . A regular sequence of 
the Joule heating is represented by a parabola-like curve, 
but the interfacial resistance is a factor of temperature 
behaviour to result in a small kink in the temperature 
distribution at the interface which is interpreted as bulk 
heating to be dominant in comparison with the interfacial 
one. 


IV. RESULTS 

In general, interfacial resistances R 1/2 may vary from 
close to zero for good metallic contacts to that of a very 
thin (Ohmic) tunnel barrier. A highly resistive inter¬ 
face can, e.g., be caused by a sample fabrication process 
in which the vacuum is broken, leading to organic de¬ 
posits. We simulate resistive E|N or N|Ng interfaces by 
modulating 2 from zero resistance to a large value. A 
large resistance of either interface turns out to enhance 
the cooling effect as long as the interfacial Joule heating 
does not dominate, i.e., when the current bias is not too 
large. 


A. Temperature profiles in a F|N|Ng pillar 


We are interested in the temperature profile in a pillar 
with equal temperatures of the two external reservoirs 
Tl — Tr = Tq, noting that the model can be easily 
extended to calculate the thermopower due to a global 
temperature difference over the device. We start with 
Tni = = Tn2 = T'j^2 = Tq as initial conditions (see 

Fig. 2]), which is substituted into Eqs. dsni) and dSlD 
to obtain the first iteration. The temperature profiles 
converge after several iterations. 

Results for F|N|Ng nanopillars are shown in Fig. [S]for 
different current densities, with temperature Tq in the 
reservoirs maintained at 300K, using parameters from 
Tables U - M for bulk and interfaces, for the case of 
all Joule heating occurring in the reservoirs. The top 
panel of Fig. [5] is for clean interfaces with AcRi = 
0.915 fllm^^ and Aci ?2 = 0.34fnm^jH. Values of in¬ 
terfacial electrical resistance are well-known parameters, 
while those of Kapitza heat conductance are not, spe¬ 
cially for F\N interfaces, and value of Kapitza heat con¬ 
ductance in latter is assumed not to be such a good 
heat conductor as compared with the second interface. 
The Joule heating is generated mainly by the relatively 
resistive ferromagnet, while the cooling takes place at 
the F|N interface, giving rise to a complex temperature 
and heat current distribution. The dotted lines for each 
curve show the average temperature in the different lay¬ 
ers Txavg govern the resistance change of the pil¬ 
lar. The bottom panel of Fig. [5] shows the tempera¬ 
ture profile in the presence of a dirty interface N|Ng 
with a 100 times larger resistance Aci ?2 = Sdlflm^. 
Gth,2 = 5.9 ■ lO^W/m^K is assumed to be reduced 
by the same ratio, while other parameters are kept the 
same. The dissipation at the dirty interface N|Ng lo¬ 
cally increases the temperature in the normal metals. 
A marked discontinuity of the temperature at N|Ng in¬ 
terface develops due to the small thermal conductance 
Gth,2- The temperature on the F-side drops from ap¬ 
proximately 298.3 K for a clean N|Ng interface to 297.7 
K for the dirty one (see Fig. 0). The increased interface 
resistance forms a barrier for the heat flow from the heat 
sinks towards the interface, allowing the region close to 
the interface to cool down more efficiently, thereby en¬ 
hancing the effective Peltier effect. 


B. Peltier cooling, Joule heating, and R-I 
characteristics 


According to Eq. (1501) the temperature profile T{x) 
is directly related to the observable resistance change. 
We compute a specific temperature profile for a given 
current bias as sketched below, which can be used to 
obtain the total resistance as a function of current that 
may be compared with experimental results. To this end 



Material 

K 

S Pf = p'f 


Co 2 MnSi 

15 

-20 0.71 

6.07 X 10""'" • {Lf/A^) 

Au (N) 

318 

1.83 

8.14 X 10"^^ • (T/Ac) 

Cu (Ns) 

401 

1.94 

6.84 X 10"“ • {Lb/AT) 


TABLE II. Thermoelectric parameters of the bulk metal lay¬ 
ers in the F|N|Ng nanopillars at 300 K: Thermal conductiv¬ 
ities K (W/mK),— and Seebeck coefficients S (fiY/K). 
Pf is the polarization of the conductivity for the ferromag- 
net while P'p is the polarization of its energy derivative. For 
lack of sufficient data we take P'p — Pf = 0, thereby dis¬ 
regarding much of the spin-dependence of the heat diffusion 
equations. Rx {Lx) are the resistances in H when thicknesses 
of the metal layers are in m. 


Material 

Gth 

S 

P 


CMS|Au 

Au|Cu 

1.8 X 10* 

5.9 X 10® 

-4 

3.5 

0.77 

0 

{dRp/dT + dRF/dT)/2 
{ORn/OT + dRNB/dT)/2 


TABLE III. Interfacial thermoelectric parameters of the 
F|N|Ng nanopillars at 300 K:'*- Interface Kapitza thermal con¬ 
ductances Gth (W/m^K) including the phonon contribution. 
S {fiV /K) is the interfacial Seebeck coefficient and P the spin 
polarization of the interface conductance. 

we linearize Eq. as: 

f) R 

^Rx^-^{Txxva-T,). (40) 

The total resistance differential is governed by the 
temperature dependence of the layer and interface resis¬ 
tances. Each bulk material layer has a specific dRx/dT, 
while the calculations establish average temperatures 
Txavg for the sections F, N and Nb respectively, as 
shown in Fig. [5] marked by dotted lines. Highly resis¬ 
tive interfaces may might affect or even dominate the 
global resistance change when Ri{ 2 ) and 9i?i(2)/9T are 
large. Our calculations include the temperatures at in¬ 
terfaces Ti( 2 ) as expressed in Eq. (1551) . For the temper¬ 
ature dependence of the bulk resistivities we adopt the 
values listed in Table |lll For resistive interfaces we aver¬ 
age dRx/dT of the two materials; This is expressed in 
N|Ng interface as represented in Table Hill 

9i?2 _ 1 / dRx dRNB\ 

dT ~ 2 \ dT ^ dT J ' ^ ’ 

while we disregard the temperature dependence of the 
resistance for good interfaces. 

In Fig. |6l the effect of inserting a highly resistive N|Ng 
interface on the R-I curves is shown for the scenario when 
the interface Joule heating is very non-local, i.e. use Eq. 
dMl)- The (effective) Peltier cooling (blue line, bottom) is 
visibly enhanced. The change in the total resistance can 
be understood in terms of the temperature distribution 
along the pillar as shown in Fig. [S] The increased interfa- 
cial resistance R 2 improves the effective Peltier coefficient 


T[K] 




FIG. 5. (Color online) Temperature distribution in a 
CMS[40nm]|Au[10nm]|Cu[100nm] nanopillar under current 
bias for the thermoelectric parameters from Tables U-EIl 
(Top) Clean F|N interfaces with resistance area of A^Ri = 
0.915 fflm^ . The brown, purple and blue lines correspond 
to applied current densities of 1,2,3 times 10^^ A/m^, re¬ 
spectively. The dotted lines for each curve show the average 
temperature in the different layers that govern the resistance 
change of the pillar. (Bottom) Temperature distribution in 
the presence of a dirty N|Ng interface with 100 times the 
electric AcR 2 = and heat resistance but otherwise 

the same parameters as the clean interface. To has a constant 
value of 300K. 


from n = 11.2 mV for a clean interface to H = 23.9 mV 
in the case of a dirty interface. Additionally, a change in 
the Peltier coefficient from H = 23.9 mV to H = 24 mV is 
reached when Eq. (1411) is implemented into this compu¬ 
tation. We should note that while the effective Peltier co¬ 
efficient is enhanced by a highly resistive interface under 
a constant current bias, it becomes a more efficient sys¬ 
tem, viz. the nanopillar requires a lower applied voltage 
in combination with more cooling effect simultaneously. 


C. Trilayer nanopillar model 

We now valuate the thermoelectric performance as 
a function of structural and material parameters of 
the nanopillars. Matching Bosu et aL’si^ samples, we 
adopt bulk (Drude) thermopowers of the leads as Sp = 
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Sc MS = — 20/xV/K for the ferromagnetic Heusler alloy 
(Co 2 MnSiTa,„„= 500 /f), Sau = 1.83^V/K for the normal 
metal N and Scu = 1.94/rV/K in normal metal Nb- 
Our model is scale-invariant with respect to the pillar 
diameter, so we cannot explain the enhanced effective 
Peltier cooling found in the narrowest pillars by the ex¬ 
periment in terms of an intrinsic size effect. However, 
smaller structures can be more susceptible to the effects 
of e.g. incomplete removal of resist material used dur¬ 
ing nanofabrication. We have disucssed above that such 
extrinsic effects do affect the thermoelectric properties 
and can be treated in our model. The interfacial ther¬ 
mopower Scms\Au and its spin polarization Ps are ba¬ 
sically unknown parameters that may contribute impor¬ 
tantly to the cooling effect in nanostructures, as reflected 
in the enhancement of the global effective Peltier coeffi¬ 
cient n = 11.2 mV for = —4/iV/K to 23.2 mV for 
Si = Scms\Au = — 30/rV/K; this case is especially rele¬ 
vant in the presence of a resistive N|Ng interface. 

The effects of an enhanced interface resistance Hci?i( 2 ) 
on the Peltier cooling can also be tested by varying it 
from that of a good intermetallic to a value correspond¬ 
ing to a thin tunnel barrier. The interface resistance 
turns out to improvement of H as long as the addi¬ 
tional Joule heating does not dominate, as illustrated 
in Fig. m Furthermore, in Fig. [3 it is plotted the 
temperature profile distribution when Joule heating is 
generated in the interfaces, setting the nanopillar with 
the same parameters of Fig. [Sj except for the mod¬ 
ulus of the electrical tunnel junction nor the Kapitza 
thermal conductance, in which both have the same ra- 


R[n] 



FIG. 6. (Color online) Resistance-current relation of 
(7445(40 nm]|T«[10 nm]|C'u[100 nm] pillars including inter¬ 
face resistances as listed in Table HIl The effective Peltier co¬ 
efficient is n = Rolp = 1.25376 Q ■ 8.95 mA = 11.2 mV (black 
line, top). The Peltier coefficient is increased to 11 = Rolp = 
2.75645 • 8.7 mA = 24 mV (blue line, bottom) when a re¬ 

sistive N|Ng interface of AcR 2 = 34fflm^ is inserted. For 
reference, the bulk Peltier coefficient is P-cms\Au = 6 mV. 


T[K] 



FIG. 7. (Golor online) Temperature distribution in the pres¬ 
ence of a dirty N|Ng interface with 10 times the electric 
AcR 2 = 3.4fflm^ and heat resistance but otherwise the same 
parameters as the clean interface. To has a constant value 
of 300K, setting the nanopillar with the same parameters as 
Fig. O but for the case of Joule heating generated at inter¬ 
faces instead. 

tio of change. It can be compared clearly a decrement 
in the performance of the cooling device for this case, 
since the Joule heating produced at the interface coun¬ 
teracts the cooling of Peltier effect. As discussed above, 
the interface resistance hinders the flow of heat cur¬ 
rent from the heat baths towards the cooling interface. 
For an interfacial resistance of AcRi = 0.915 fHm^ and 
Aci ?2 = 0.34fHm^, the total Peltier coefficient reaches 
a value of ^cms\Au\Cu = 11.2 mV, matching parame¬ 
ters from Tables Hi] and IIIIl where this result from this 
theoretical model is close to experimental ones. A lin¬ 
ear dependence of the Peltier coefficient was found when 
varying the interface resistance area AcRi from 0.915, 
9.15 and 91.5 (fHm^), resulting in Peltier coefficients 
ncMS|Au|Cu of 11.2, 13.49, and 31.61 mV, respectively. 
By contrast,when the interfaces are clean and Joule heat¬ 
ing is suppressed (assuming -I- = oo), the Peltier 

coefficients increase to 11.28, 14 and 42mV for the same 
interface resistances for the best case when Joule heating 
is all produced in the reservoirs. 

Since our calculations take the spin degree of freedom 
into account the spin accumulations and spin currents 
along the nanopillar are byproducts of the calculations. 
In contrast to |Pf| < 1, the spin polarization of the 
derivative of the conductivity —oo < Pp < oo. When 
Pp < Pp the spin contribution to the cooling power 
is proportional to the spin accumulations as expressed 
in Eqs. ((5511 and Eq. ((571) . A Peltier coefficient of 
n of 11.2 mV with parameters from Tables HD and HID 
is increased by a factor 2 when Pp = —20. However, 
if Pp > Pp, the spin degree actually generates heating 
thereby reducing the cooling power. 

We also studied the dependence of the effective cooling 
on layer thicknesses Lp,L and Lq. The Joule heating 
dominates for a critical current bias Ic that decreases 
with increasing Lp. When the thickness of F=CMS is 

















10 


reduced from 40 nm to approximately 5 nm, If improves 
slightly from 11.2 to 12.5 mV. The optimal thickness of 
the ferromagnetic him is Lp ~ \p. The normal metals 
do not signihcantly contribute to the cooling since their 
Peltier coefficients are relatively small. 

Finally, slight enhancements of the Peltier coefficient 
could be achieved by including in the analysis an external 
heat current JqCxt, which is depicted in the left hand of 
Fig.H which forms part of an extension of the nanopillar 
that could lead the head current towards a further reser¬ 
voir so that Tp < To, to result in a slight enhancement 
of the Peltier effect. This makes a more sophisticated 
model, but we leave it for a future study. 


V. SUMMARY AND CONCLUSIONS 


This paper is motivated by the observed enhancement 
of the cooling power in magnetic pillars when the cross 
section was reduced to the nanoscale.— We develop a 
realistic spin, charge, and heat diffusion model to inves¬ 
tigate the roles of spin-dependent bulk and interface scat¬ 
tering contributions. We analyzed the (apparent) cool¬ 
ing power and the conditions to maximize the effective 
Peltier effect. 

We demonstrate that very thin (Ohmic) tunnel junc¬ 
tions can improve the cooling power of devices as appar¬ 
ent in the shift of R{I) parabolas. On the other hand, 
the spin degree of freedom that was thought to be essen¬ 
tial in CMS materials appears to be less important for 
conservatively chosen parameters. However, the mate¬ 
rial dependence of key parameters is basically unknown. 
The parameter Pp,i.e. the spin polarization of the spec¬ 
tral asymmetry of the conductance, turns out to play an 
important role. This parameter may become arbitrarily 


large when d 





= 0 or, for interfaces 


d (g; + G^) 



= 0, which does not seem to be an 


exotic condition and we recommend a systematic search 
for such materials or material combinations. Our results 
also indicate that interfacial parameters such as the in¬ 
terface Seebeck coefficients >S'i( 2 ) play a very signihcant 
role in the thermoelectric characteristics of multilayers 
and may not be disregarded when validating their perfor¬ 
mance and in agreement with Hu et al^ affirming that 
the value of this coefficient is even larger that the con¬ 
ventional one, that in the present model both contribute 
simultaneously in the cooling effect. 

While the experiments up now have been analyzed in 
a simplistic model for the compensation current at which 
heating and cooling effects cancel, we established a dis¬ 
tributed model of currents and temperatures. The com¬ 
puted temperature prohles along the nanopillar estab¬ 
lished that the cooling is not homogeneous, but heating 
and cooling coexists in different locations of the sam¬ 
ple. The current-dependent resistance only a very crude 
thermometer that is not a reliable measure for a cooling 
power that could be of practical use. 


We hnd that it is possible to selectively cool a ferro- 
magnet by a few degrees simply by a current flow in the 
right direction. This could be an important design pa¬ 
rameter for STT-MRAMs. The writing of a bit of infor¬ 
mation by a switching event of the free layer in a memory 
elements is accompanied by signihcant Joule and Gilbert 
heating. Applying a small bias current after the magne¬ 
tization reversal can assist a quick return to the ambient 
temperature. 

Our model is scale invariant with respect to the pillar 
diameter and does not provide and intrinsic mechanism 
for the observed size dependence of the Peltier effect. In 
principle, extrinsic effects should exist. The large fluctu¬ 
ations observed in the experimental results indicates sig¬ 
nihcant disorder in the smallest nanopillars. One source 
of the problems can be the need to break the vacuum 
during sample fabrication. The effect of pollutants at an 
interface are then likely to be more serious for smaller 
pillars. We found indeed that by modelling interface as a 
thin tunnel junction enhances the apparent Peltier coef- 
hcients by suppressing the heat currents howing into the 
pillar from the reservoirs. However, the record cooling ef¬ 
fects observed for some of the narrowest pillars appear to 
be beyond the effects that can credibly be modelled, and 
we cannot exclude the possibility that something more 
interesting is going on. 

Several effects are beyond the present model ap¬ 
proach. Size quantization is not expected to be important 
in metallic structures at room temperature, but could 
play a role in heterogeneous materials disordered on a 
nanometer-scale. Spin waves and magnons, i.e. exci¬ 
tations of the magnetic order parameter, affect thermo¬ 
electric properties. The magnon-drag effect^ enhances 
the Seebeck coefficient . The longitudinal spin Peltier 
effect^ discovered for bilayers with magnetic insulators 
should also exist in metallic structures: the spin accu¬ 
mulation in the normal metal generates a heat current 
that comes on top of the heat currents discussed here. 
It is not clear, however, how and why these effects be¬ 
come so strongly enhanced in the nanopillars addressed 
experimentally. More experiments on even smaller and 
more reproducibly fabricated nanopillars, preferably fab¬ 
ricated without breaking the vacuum, are necessary in 
order to provide hints on what is going on. 

We conclude that the Peltier effect in magnetic 
nanopihars with diameters > 100 nm appears to be well 
understood, but that the enhanced values for narrower 
ones are to date only partly explained. In order to em¬ 
ploy the large observed effects, more experiments are nec¬ 
essary in order to shed light on the underlying physical 
mechanisms. 
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